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plicability of the replica approach to description of spectral fluctuations in the context of random 
matrix theory and beyond. The present paper, concentrating on invariant non-Gaussian random 
matrix ensembles with orthogonal, unitary and symplectic symmetries, aims to demonstrate that 
both the bosonic and the fermionic replicas are capable of reproducing nonperturbative fluctuation 
formulas for spectral correlation functions in entire energy scale, including the self-correlation of 
energy levels, provided no a-model mapping is used. 
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I. INTRODUCTION 

Energy level statistics of electrons in a random poten- 
tial is one of the basic issues in the theory of disordered 
conductors. It is also a classic example of the problem 
whose proper statistical description calls for a use of non- 
perturbative methods of averaging over disorder realiza- 
tions. The latter procedure is not trivial because the 
spectral observables depend in highly nonlinear fashion 
on a random Hamiltonian, thus making the calculation 
of the ensemble averages very difficult. Within the single 
electron picture, it is necessary to average a product of 
resolvents, (r|(z — 7i)~^ |r') , that admits an exact rep- 
resentation through a ratio of two functional integrals 
over auxiliary fields which may consist of either com- 
muting [bosonic] or anticommuting [fermionic] entries. 
This representation, however, is still not convenient for 
nonperturbative averaging due to an awkwardrdenomina- 
tor. To get rid of it, Edwards and AndersonEJ proposed 
a replica method that substitutes the initial disordered 
system by the n identical nonintcracting systems. Such a 
trick removes a randomness from the denominator and, 
therefore, makes it possible to perform a nonperturba- 



tive averaging over the random potential from the very 
beginning for all integer n which must be set to zero 
at the appropriate stage [by the analytic continuation 
n — > 0]. Depending on the origin of the auxiliary fields, 
the resulting, effective field -theory is defined on cither a 
noncompactEl or a compactB manifold. j-. 

An alternative approach was invented by EfetovQ who 
introduced an auxiliary field with an equal number of 
bosonic and fermionic entries that obviates the need 
for replicas. This supersymmetry formalism turned out 
to be very fruitful in the theory of disordered conduc- 
tors and beyondd. Along with othec. fundamental re- 
sults, it has led to establishing a linkau between the the- 
ory of disordered metals and the random matrix theory 
[RMT] founded by WignerQ, thus confirming on micro- 
scopic grounds a much earlier conjecture by Gorkov and 
EliashbergEl. To be precise, it was demonstratedo that the 
'density-density' correlation function Rp^s) of disordered 
conductor in the ergodic limit is solely determined by 
the global symmetries of the system and coincides with 
the RMT predictions. For the three symmetry classes 
- orthogonal {/? — 1), unitary {(3 — 2), and symplectic 
{(3 = 4)- the 'density-density' correlation functions are 
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s being the distance between the eigenlevels measured in 
units of the mean level spacing A. The range of validity 
of Eqs. (LI) is s ^ g, where g » 1 is the dimensionless 
Thoulcss conductance. 

While the success! of the supersymmetry approach 
has beeipiiyery impressive, the alternative replicated field 
theoriesH'U [which, besides, are suitable for a treatment of 



disordered systenjs in presence of interactions^] seemed 
to be restrictedll3 to p^producing the perturbative dia- 
grammatic expansionsEiJ only. A detailed analysis of this 
problem, has been given by Verbaarschot and Zirnbauer 
in RefJl3, where the failure of the replica method [in its 
(T-model formulation] to correctly account for all nonper- 
turbative contributions has been attributed to the fact 
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that the replica trick is mathematically ill founded, and 
the knowledge of the replicated partition function for all 
integer n need not be sufficient for extrapolation to the 
limit n ^ Q which, generically, can also be nonuniqueEj. 

Surprisingly enough, in the very recent publicationQ, 
Kamenev and Mezard have formulated a variant of the 
fermionic replica method that allowed them to rederive 
Eq. dl.lbD , s » 1, in the context of RMT fa^the Gaus- 
sian Unitary Ensemble [GUE] . The approachliJ has relied 
on the combination of the exact integration over the an- 
gular degrees of freedom of the replicated cr-model, by us- 
ing the Itzykson-Zuber integraltj, and the approximate 
saddle point calculations in the n-dimensional space of 
remained eigenvalues A — diag(Ai, . . . , A ^). In order to 



correctly reproduce the result Eq. (1.1b) for GUE in the 
large-s limit, one is forced to take into account the non- 
trivial saddle points which break the replica symmetry 
between the different eigenvalues A;, 1 < i < n. On the 
formal level, this implies that the n-component vector of 
eigenvalues ascribes a nontrivial internal structure which 
must be kept in the limit n ^ of vanishing number of 
integration variables! Very similar structure of the vec- 
tor saddle Ppi|Pt] manifolds is present in two subseowent 
publicationstSim that extended the results of Ref.ll3 to 
two other [orthogonal a»d symplectic] syipHietry classes 
both in the RMT limitO and beyond ittSlla. A-m 



cal discussion of the procedures employed in Refs.liJIl^ 
has been given in Ref.t2l where it has been shown that 
there e xists .. so me additional input in the calculational 
schemestSEJO which is equivalent to introducing the 
causality information. A conjecture was nut , f prwardEj 
that the fermionic replicas in elaborationlljEJO are per- 
turbatively equivalent to the supersymmetry method- 
with being iSr^mall parameter. It was also explainedE2l 
that the resultal^O are 'semiclassically exact' for the uni- 
tary syntmetry class owing to the Duistermaat-Heckman 
theoremE3. 

There is one common feature inherent to the above 
mentioned replicated calculations: All of them, start- 
ing with either bosonic or fermionic replicas, reduce the 
problem of evaluating the spectral density-density cor- 
relation function to the effective cr-model defined in the 
configurational space of replicas only. Correspondingly, 
the 'coordinate' space disappears from the theory, and 
the memory of the initial [matrix] Hamiltonian is only 
retained in the symmetry of the auxiliary Q-fields. An 
immediate consequence of such a treatment is that the in- 
tegral representations for the spectral observables appear 
to involve the integration measures which become sense- 
less as n ^ 0. This can be seen, for instance, from the 
Verbaarschot-Zirnbauetintegral for the two-level corre- 
lation function derivedtHI from the fermionic replicas for 
the GUE in the large- hmit: 
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[Here, r — tts in notations of Eq. (1.1b)]. Obviously, one 
cannot simply take the limit n — > neither in the above 
integral representation nor in the result of the integration 
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because of the vanishing dimensionality of the matrix un- 
der the sign of the determinant. Therefore, one sh ould 
know how to evaluate the integrals of the type Eq. (1.2) 
in the closed form for arbitrary integer n in order to an- 
alytically continue the resulting expression from integer 
n and perform the replicaJimit n — + afterward. One 
such way proposed in Ref.llZl reduces Eq. (1.2), in the do- 



main 7' ^ 1, to the expression which coincides identically 
with the one obtained by means of the 'replica symmetry 
breaking' [RSB]. This coincidence signals that the RSB 
may be viewed as a useful language for approximate com- 
putation of the integrals defined on the problematic [in 
the above mentioned sense] measures. In the course of 
exact calculations, the notion of RSB must not appear 

The studiesEi'Q'Ei'liBEi call for further clarification of 
the question that has a conceptual importance and con- 
cerns the applicability of the replica method to a non- 
perturbative description of the eigenlevel correlations in 
general, and those on the short distance scale, in partic- 
ular. The central issue to be learnt here is whether the 
replica method itself faces some internal insurmountable 
obstacles, or existing difficulties are to be attributed to 
the particular computational schemes. We believe that 
in order to answer this question, one has to depart from 
conventional routes toward exact calculational schemes 
within the replica framework. 

While it is unknown so far how the exact replica treat- 
ment could be realized for realistic disordered systems, 
the present paper aims to demonstrate that there is one 
particular but conceptually important situation where 
the exact computation is possible for both the fermionic 
and the bosonic replicas.r-.This is the case of invariant 
random matrix ensemblesu with a general non-Gaussian 
measure 
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specified by the confinement potential ^[H] and the 
Dyson index /3 = 1,2 and 4 that labels the symmetry 
of the matrix Hamiltonian H. 

A distinctive feature of our treatment that must be 
emphasized from the very beginning is that we do not 
use a cr-model representation for the replicated parti- 
tion function averaged over the ensemble Eq. (1.4) of 



random matrices. Instead, a different technique is used 
for exact averaging of the product statistics that com- 
prises the replicated partition function as a particular 
case. This technique generates especially useful repre- 
sentations for the spectral observables, thus allowing one 
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to avoid the problem of the ill defined integration mea- 
sure. As a direct consequence of that, we construct [in 
a simple and straightforward way] a well defined replica 
limit n — > in expressions for the spectral observ able s 
that turn out to identically coincide with Eqs. (1.1). 
This coincidence, observed for arbitrary s, is the main 
outcome of the study. 

The following program is carried out: 

(i) First, the replicated partition hmction of the form 
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X det" [(^' - H)2 + 7'2] (1.5) 

is introduced in accordance with the very concept of the 
replica method that boils down to constructing such a 
generating function [e. g., for the 'density-density' cor- 
relator] which does not contain a random denominator 
at the expense of the artificial appearance of replica pa- 
rameter n in Eq. (1.5). 

(ii) Second, the replicate d pa rtition function is aver- 
aged over the ensemble Eq. (|l . 4|) .of xapdoaxjuatrices . At 
this stage, we depart from Rcfs£3'oE30'E3 and do ex- 
act averaging without appealing to a nonlinear cr-modcl. 
This is the most crucial step in our derivation: It allows 
us to retain control over the analyticity of the resulting 
expressions considered as a function of n, the number of 
replicas. 

(iii) Third, the 'density-density' correlation functicp,is 
recovered by means of the double limiting procedureEa 

{pn{E)pn{E')) ^ •■ ^ 



lim lim , ^„ 
7,7'^On^±o (27rn)"^ 

X d^dy{Zl^,{E,E')) 



(1.6) 

and found to follow Eqs. (1.1) in the spectrum bulk after 
appropriate rescaling of spectral variables. [The limits 
n — > ±0 correspond to the fermionic and the bosonic 
replicas, respectively.] 

The paper is organized as follows. In the auxiliary Sec- 
tion II, the notion of the product statistics is introduced, 
and a technique to average the product statistics is pre- 
sented. With minor modifications, this sectioiL. closely 
follows the recent paper by Tracy and WidomEj whose 
impact we acknowledge. The results of Section II, repre- 
sented in the form of the Grassmann functional integrals, 
are us ed i n Section III where the fluctuation formulas 
Eqs. (LI) are obtained. Conclusions are presented in 
Section IV. 



II. PRODUCT STATISTICS IN RANDOM 
MATRIX THEORY 

Let us introduce the notion of the product statistics F 
defined on the eigenvalues {A} of N x N random matrix 



H of appropriate symmetry: 
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Here, /(A) is an arbitrary function of the variable A such 
that the integral J dXX' f[X)e~^^^^^ exists for positive 
integers j. 

This section aims to establish exact formulas for the 
average of the product statistics 
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over the distribution function Eq. ( |1.4| ) for all three sym- 
metry classes /3 = 1, 2 and 4. We notice that along with 
utility of these formulas for computing the average of the 
replicated partition function, they may find applications 
in evaluating the average of other product statistics-ep-, 
countered in description of a number of observablesEllE^ 
in mesoscopic physics. 



A. Unitary symmetry: (5 = 2 



Due to invariance of the distribution function P/v[H] 
aX (3 — 2 under the unitary transformation U(iV), the 
average {F) can be written down as TV-fold integral 
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where the proper normalization is to be restored 
on the later stage. Here, the diagonal matrix 
A = diag(Ai, . . . , Aat), the integration measure I?A = 
]\m=i '^^rn, and AAr(A) = ]\i<i^j<N{K~ Xj) is the Van- 
dermonde determinant. 

Integration in Eq. ( ^.3[ ) can be performed if we 
notice that the Vandermonde determinant AAr(A) = 
det[Af~^] oc det[Pj_i(Ai)], where Pj{X) is a set of polyno- 
mials that we choose to be orthonormal with respect to 
the measure e~^^^'^^(iA. In terms of the 'wave functions' 
VjW — e^^^'^'-'Pj (A) this orthonormality is expressed as 



dXipi{X)(p.j{X) = 5i. 



(2.4) 



The integral Eq. (2.3) can be brought to the form 
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Y[_ dXmf{Xm) det[(^j_i(Aj)]i<ij<Ar det[(^j_i(Ai)]i<ij<Ar. 



(2.5) 
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Integrating the product of two determinants in accordance with the formulae^ 



/ JJ dX,n det[uj{\i)]i<i,j<N det[wj{\i)]i<i,j<N = iV! det ( / d\u.i,{X)Wj{X) 1 

\m=l J \J / l<t,j<N 



(2.6) 



that is proven via the Laplace expansion of the two de- 
terminants, we reduce Eq. (|2.5|) to 



(F) =dct( f dA/(A)(p,_i (A)(p,_i (A)) (2.7) 



that enjoys a correct normahzation owing to Eq. (2.4) 
and constitutes the first preliminary result to be used in 
the following Section. 



B. Symplectic symmetry: (3 — 4: 

In the case of the symplectic symmetry, Sp(iV), the 
average is given by 



(2., 



where the correct normalization is to be found yet. De- 
spite of appearance of |AAr(A)|^, that reflects the self- 
dual quaternion nature of the matrix H for /3 = 4 in Eq. 



(1.4), the exact integration in Eq,— (2.8) is still possible. 
First, we make use of the identitynZI 



|AAr(A)|4 = det 



(2.9) 



det [P,_i(AO,P;_i(A,)], 



<i<N,l<j<2N 



to rewrite Eq. (2.8) as 
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In both equations above, PjiX) is a set of polynomials 
which we choose to be orthonormal with respect to the 
measure e~^^*^'^^dA. Second, thcjatcgrati on in Eq. ( 2.10| ) 
can be performed via the analogO of Eq. ( |2.6| ) that reads 



J i n dKnj dct[Uj(AO,U>j(Aj)]^<^<^_^ 



<j<2N 



^{2Ny.pf[ I dX[u,{X)wjiX)-u,iX)w,iX)]] . (2.11) 

l<iJ<2N 



Here, the notation 'pf stands for the pfaflianE3 [its square is determinant]. The identity Eq. ( 2.11 ) enables us to 
simphfy Eq. (2.10) to 



pf 



/ dA/(A)[^._i(A)v:>;_i(A) - ^,_i(A)^^_i(A)]) 

/ l<i,j<2N 



(2.12) 



where irrelevant overall constant has been omitted, and 
the 'wave functions' fj{X) = e^^^^^^ Pj{X) o beyin g the 
orthonormality relation in the form of Eq. (2^) have 
been introduced. Combining the above results together 



and properly restoring the normalization, we derive the 
following exact pfaffian representation for the averaged 
product statistics: 



pf (/ rfA/(A) (A)(^;_i (A) - ^j-i{XM_i (A)]) , 
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Equation ( 2.13 ) is the second preliminary result needed 
in what follows. 



C. Orthogonal symmetry: /3 = 1 

For (3 = 1, the pfafRan representation for the averaged 
product statistics similar to Eq. (2.13) can also be ob- 
tained. To avoid unnecessary complications, we choose 



the real symmetric matrix H to be of an even dimension, 
2N X 2N , so th at t he upper index in the products in 
Eqs. ( |2T| ) and ( pl| ) is equal to 27V. Then, due to the 
0(27V) invariance of the distribution function P27v[H], 
we encounter the following 2Af-fold integral 



(P)oc / I?Ae-"-^(^)|A2Ar(A)|det[/(A)]. (2.14) 



It is useful to reduce Eq. (2.14) to the formi 
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X n (^J-^O, (2.15) 

l<i<i<2N 

where the integration runs over the ordered domain 
Ai < ... < A2Ar, and therefore the modulus of the 



Vandermonde determinant has been relaxed. Recalling 
that the Vandermonde determinant is proportional to 
det[Pj_i(A,,)] with 1 < i,j < 2N and choosing Pj{X) 
to be the polynomial orthonormal with respect to the 
measure e' we are able to pieijform integration in 

Eq. ( 2.15 ) by means of the identityEa 



J '^^"^ det[iij(Ai)]i<jj<2Ar = pf ^A J dX'sgn{X' - X)ui{X)uj{X')j 



<i.j<2N 



where sgn(A) = A/|A|. The integration results in 



pf( fdX f dX'e{X - A')/(A)/(A')¥'»-i(A)^,_i(A')) 

\J J / l<i,j<2N 



where irrelevant overall constant has been suppressed, 
and the 'wave functions' (pjiX) = e'^'^^^^'^ Pj{X) obey- 
ing the orthonormality relation in the form of Eq. (2.4) 
havei-appjsared. Also, to meet the standard notations of 
RefsEl'E2l, we have introduced 



e(A) = -sgn(A). 



(2.16) 
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Restoring the proper normalization, one arrives at the 
third preliminary result: 
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pf (/ / dXdX't{X - A')^.-i(A)^,-i (A')) 



(2.19) 
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Notice that Eqs. (^, ( |2.13| ) and ( |2.19| ) are exact and 
refer to an arbitrary product statistics F that does not 
affect the convergence of the integrals under the signs of 
'det' and 'pf. 



III. EIGENLEVEL CORRELATIONS FROM THE 
REPLICA METHOD 

A. Unitary symmetry: /? = 2 

Exact determinantal expression for the averaged repli- 
cated partition function {Z^ ^,{E, E')) in the case of the 



unitary symmetry is readily obtained from Eq. (2.7) 
upon specifying 



f{X)^[{E-Xr+j']-[iE'-Xr+^'']\ (3.1) 



see Eq. (IJ). It should Ja£Ljppticed that, contrary to the 
tJ-mo del representatiorOo'EZlEj [see, also, Eqs. (1.2) 



and (1.3)], our determinantal form is especially useful 
to perform the replica limit n — > ±0, with -f and — 
corresponding to the fermionic and the bosonic replicas, 
respectively. However, Eq. (2/7) as it stands, is not con- 
venient for the large-TV analysis of the 'density-den sity ' 
correlation function. For this reason, we put Eq. (2.7) 
into equivalent form of the functional integral over the 
scalar Grassmann field 
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with the measure 
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(3.3) 
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where {xk} are N Grassmann complex variables_and t/j 
denotes the complex conjugate. One can verifyEB that 
the choice Eq. ( |3.2| ) enables us to rewrite Eq. (2.7) in 
the form of the functional integral 



{F}^ D 



that solely involves the scalar kernel! 

N-l 



i 



KN{x,y) = ^ Lpk{x)ipk{y) = Kiqijj.x) 

k=0 

obeying the integral relationship 

dzK^ix, z)Kn{z, y) = Kn{x, y). 



(3.5) 



(3.6) 



In turn, the universal properties of the scaiaip-kernel 
Kn{x, y) in the large- A'' limit are well studiedtilH. 
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Substituting Eq. (3.1) into Eq. (3.4), differentiating 
the latter with respect to 7 and 7' and computing the 
limits n — > ±0 and 7, 7' one gets access to the 
'dens ity- density' correlation function in accordance with 
Eq. (1.6). Since we are interested in the correlations in 
the energy range where E — E' may take the values of 
order of the mean level spacing A^v, we introduce the 
scaled variables S — E/An and 5" = E'/A^, and re- 
strict ourselves to the vicinity of the center of the eigen- 



value support, where S, S' ^ 0{N^). In the large- iV 
limit, the rescaled Eq. ( [Lq ) is then reduced to 



lim lim 



1 



X d^dy{ZlyiSAN,S'AN)). 



(3.7) 



Straightforward calculations base d on Eqs. (3.7) and 
~4|) with /(A) given by Eq. (O) yield 



{p{S)p{S')) =diS ~S')A% J daKNiaAN,SAN){maANMS^N))) 

+ ^%J J dada'KN{aAN,SAN){{^{ crAN)i^iSAN)ip{(T'AN)ip{S'AN)))KN{cr'AN, S'An), (3.8) 
where the double angular brackets ((. . .)) stand for the Gaussian average 



(3.10a) 



In Eq. (p]q), {{tpip)) a-nd {{ipipipil^)) may be computed via the Wick theore: 



((■0(0-1 Ajv)V'(o-2Ajv))) = KNicTiAN, (T2An), 
{{i>{alAN)'ip{a2AN)^Pia■3AN)i>ia■4AN))) = KNi<JiAN,(J2AN)KN{cr3^N,cr4AN) 

~ KNi<JiAN,a^AN)KNicT3AN, ctsAat). (3.10b) 



As far as we are concerned with the spectrum bulk, 
one has AnKnISAn, S'An) = K{s) with 



K{.s) = 



sin(7rs) 



ITS 



(3.11) 



s = IS" - 5" I, by the universality argumentf£ilH'El. This, 
equation holds for the strong confinement potentialo 
only. 



Combining now Eqs. (3.J), ( 3.10 ) and ( 3.11 ), we come 
down to 



{piS)p{S')) 



l-K'(s) 



6{s) 



(3.12) 



which is clearly identical to Eq. (l.lt). This is the cen- 
tral result of the subsection. 



It must be stressed that, by derivation, Eq. ( 3.12| ) 
is valid in the bulk of the spectrum without the restric- 
tion s ^ri] inherent to the approximate calculational 
schemesoEZl. The self-correlation of energy levels is also 
reproduced. 



B. Symplectic symmetry: /3 = 4 

For technical reasons dictated by the pfaffian structure 
of Eq. (2.13) derived for (3 — 4, it is useful to deal with 
{Z'^yiETE^}^ instead of {Z:^y{E,E')). In such a case, 
the 'density-density' correlation function is seen to be 
generated by the replica limit 



{pn{E)pn{E')) 



lim lim 



2 7,7'-*0 n^±0 (27m)2 



d,dy{Zl^,(E,E')f - -d,{Zl^,(E,E'))%,{Z-^^,{E,E')f 



(3.13) 



with (Z" {E, E')) being the ratio of two determinants. 



Eq. ( |2.13D ; the function /(A) is given by Eq. (fU|). 

To facilitat e the large-A^ analysis, we rewrite the 
squared Eq. (2.13) in the equivalent form of the ratio 
of two Grassmann functional integrals: 



{Ff = ^ L^-^ = — , . (3.14) 

j jji^ ^ J ^ dxdyi'(x)K„(x,y)^(y) 



Here, ^'(x) is the two-component Grassmann field 

2^ 2N , . 

fc=l fc=l ^ ^ ' / 



and *(a;) = '^\x)L, with 



L 



-1 

1 



(3.16) 
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Also, the integration measure is 

2N 



M,k = ydA[^,_i(A)^l._i(A)-^fc_i(A)^;_i(A)], (3.21) 



(3.17) 



fe=i 



The choice Eq. (3.15) naturally leads to appearance of 
the 2x2 matrix kernel 



Kll{x,y) K]^{x,y) 



(3.18) 



whic h is a close analog of the scalar kernel KN{x,y), 
Eq. (3.5), arising in random matrix ensembles with V{N) 
s ymm etry. In terms of the wave functions Lpk{x) of Eq. 
(2.13), the (a,/?) entry of the matrix kernel is given by 



j,k G (l,...,2iV). Equation ( 3.19 ) coincides witi the 
one found for the first time by Tracy and Widomc3. As 
is the case of U(iV) symmetry, the large:^iV^hehavior of 
the matrix kernel Kn is known in detailEifE3l23. 

The representat ion Eq. ( ^j.l4 ) is a bit less trivial as 
compared to Eq. (3.4). To verify the former, one should 
observe that (j)k{x) in Eq. ( 3.15 ) is the right eigenvector 
of the matrix kernel KN{x,y): 



(3.22) 



27V 



Kf(x,y) = (-1)^-1 ^ ^f_r\x)^,,k^t-'Hy), 

3,k=l 



(3.19) 



[a, /3 = 1, 2], so that the integral identity 
dzk^ix, z)Kn{z, y) = k^ix, y) 



dyKN{x,y)(j)k{y) = (t>k{x), 
while 4>k{.x) — 4'\{x)L is its left eigenvector: 
dx'(j)k{x)kN{x,y) = '4)k{y)i 



(3.23) 



(3.20) withfce (0,...,2Af-l). 

In what follows we are interested in the 'density- 
holds. Here, {x) denotes the derivative of a~th order density' correlation function in rescaled energy variables 
with respect to x, and ^ijk — {M~^)jk, where 2A^ x 2A^ S — E/A n and S' = E'/An, when the large- hmit of 



antisymmetric matrix M has the entries 



Eq. (3.13) transforms to 



{piS)piS')) = 



A 



N 



lim lim 



2 7,7'^On^±o {2-Kny 



1 



d^dy{ziy{SAN,s'AN)f--d^{ziy{SAN,s'AN)fdy{z:;^^,{SAN,s'AN)r 



(3.24) 



Substituting Eq. (3.14) with /(A) given by Eq. (3.1) into Eq. (3.24) leads us to the following representation: 

1 " ' 



{p{S)p{S')) = - 



5{S-S')A% [da J2 k%^'{'J^N,SAN){{^Sj"iaAN)^^iSAN))) 

2 2 

dada' J2 E K';f{'^^N,SAN)kl\a'AN,S'A, 

a, 13=1 7,(5=1 

X (((*"(crAAr)^'''(S'AAr)¥^(a'AAr)1''^(S"Ajv))) 



^{{'f''{aAN)^^{SAN))){{^\a'AN)^\S'AN))) 



Here, the double angular brackets ((. . .)) denote the Gaussian average 

J *]( )e'^" / / '''^'i'^'*('^'^N)-R'iv(<TAjv,o-'Ajv)*(cr'Ajv) 



((■••)) 



J ^Pje'^" ^ ^ d(T(i(T'*((TAjv)i?iv(<TAjv,(T'Ajv)1'(cr'Ajv) 



(3.25) 



(3.26) 



The averages of the type ((vl'"^'^)) and {{'i'"'^^'^'^ '^^)) in Eq. ( 3.25D are evaluated by means of the Wick theorem 
[see Appendix for details]; below we display the result: 



{{^^{aiANWia^A^))) = ifj^"(a2AAr, Aat), 
{(W"(aiAjv)*''(a2Ajv)W^(a3Ajv)*'(a4A^))) = i:^"((72AAr, fxi Aw)A'^^(a4Ajv, ^TjA^) 

- if^"(a4Aw,CTiAw)A:^'^(a2Ajv,a3A^) 



(3.27a) 
(3.27b) 
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Inserting Eqs. (3.27) into Eq. (3.25), and making use of Eq. (3.20), one obtains: 



{p{S)p{S')) - - [d{S - S')ANtr[kN{SAN, SA^)] 

+ A% (^tr[kN{SAN,SAN)]tT[KN{S'AN,S'AN)]~tr[KN{SAN,S'AN)KN{S'AN,SAN)] 



(3.28) 



Here, the trace 'tr' acts in the 2x2 space of the matrix 
kernel K^. 

In the large-TV limit, the kernel kN{SAN,S'AN) 
taken in the bulk of the spectrum Isi/hi(ih.is of our primary 
interest] obeys the universal law£3'E30 

AMkU{SAN,S'AN) = ANk]^(SAN,S'AN) 

= K{s), (3.29a) 

A%kj,^SAN, S'An) = K\s), (3.29b) 
k^^\SAN,S'AN) - 



dtK{t), 



provided the confinement potential is strong 
s — \S — S'\, and 

sin(27rs) 



Kis) 



2tts 



(3.29c) 
Here, 

(3.30) 



Equation ( ^.31 ), derived fr om t he replicas, is seen to iden- 
tically coincide with Eq. (1.1c) for arbitrary s. It repre- 
sents the main result of the subsection. 



C. Orthogonal symmetry: l3 — 1 



As is the case of the symplectic symmetry, it is use- 
ful to deal with (Z^ ,^,(£', ^;'))^ so that the 'density- 
density' correlation function is defined by the replica 
limit Eq. ( |3.13 ) with N being replaced by 2N, and with 
{ZJ^ ^,{E, E')Y given by the ratio of two determinants as 



Combining Eqs. ( |3^ ), ( |3^ ) and ( |33c| ), we come 
down to 

dK{s) 



{p{S)p{S'))^l-K%s) + 



ds 



dtK{t)+5{s). 



(3.31) 



follo ws fr om Eq. (2.19); the function /(A) is specified by 
Eq. (13). 

In order to simplify the large-TV analysis, we have to 
express {ZJ^ ^,{E, E'))'^ in terms of the Grassmann func- 
tional integrals. One can verify that the squared averaged 
partition function admits the representation 



(3.32) 



Here, the integration measure £'[^',5'] is given by Eq. 
(3.17). The two-component Grassmann field ^ is defined 
as 



with 



2N 



Q'U'.v) = (-1)'' Y. [^'"-Vi-il(i)fj*[^'-'V*-il(!/), 



(3.33) 



(3.36) 



with vl'(x) — '^\x)L; 2x2 matrix L is given by Eq. 
( |3.16 ). Al so, £ denotes the integral operator associated 
with Eq. ( 2.18[ ); it acts on some function g{x) in accor- 
dance with the rule 



[a,/3 = 1,2]. Here, the matrix pjk = [M ^)jk, where 
2N X 2iV antisymmetric matrix M has the entries 



dxdyipj-i{x)e{x - y)ipk-iiy), 



(3.37) 



dye{x - y)g{y) 



(3.34) 



provided the integral is convergent. 



The representation Eq. (3.32) involves the 2x2 ma- 
trices Q2N and K2N which are related to each other as 



hN[x,y) = K2N{x.,y) + 











eix-y) 



(3.35) 



j,fce (1,...,2A^). 

Let us point out that the raa-trix K2N{x,y) is noth- 
ing but the 2x2 matrix kernelHj describing energy level 
correlations in random matrix ensembles with orthogonal 
symmetry. However, contrary to Eq. (3.20), the kernel 
k2N{x, y) obeys the more complicated integral identity 



dzk2N{x,z)K2N{z,y) = K2N{x,y) 
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~]^[k2N{x,v),K], (3.38) 

where [...,...] stands for commutator, and the matrix 
A = diag(+l, —1). The auxihary matrix Q2n{x, y) satis- 
fies 



dzQ2N{x, z)Q2n{z, v) = 2Q2n{x, y). 



(3.39) 



An advantage of the representation Eq. ( p. 3^ ) is in the 
fact that the Gaussian averages of the type ((^ ^'^}) and 



arising when implementing the rephca 
limit may easily be expressed in terms of the large- 
matrix kernel K2N] the universal properties of t 
for the strong level confinementEj are well studiei 

Restricting ourselves to the spectrum bulk, the 
'density-density' correlation function in rescaled vari- 
ables S — E//S.2N and S' = E' / A2N is given b y the 
replica limit Eq. ( ^.24| ). Substituting there Eq. ( |3.32| ) 
with /(A) specified by Eq. (3.1) and taking into account 
Eq. ( 3.35 ), we obtain 



{p{S)p{S')) = ^ [e{S - S') [((*'(5'A2Ar)vI/l(5A2Ar))) - {{•^\S^2n)^\S'^2n))) 

+ 5{S-S') I dae{S-a) \{{^\a^2N)^\Sl^2N))) - {{^\s1^2n)'^\^1^2n))) 



dada'eiS' ^a)e{S-a') ((* (5'A27v)*'(aA2jv)* {S 1^2n)^\o' ^2n))) 



+ ((*'(aA2Ar)*l(5'A2Ar)^'(a'A2w)*\5A2Ar))) 

- ((*'(5'A2Jv)*'(^TA2Ar)^'(a'A2Ar)*l(5A2Ar))) 

- ((*'(aA2Ar)*l(5'A2Ar)^'(5A2Ar)*l(a'A2Ar))) 

- i (((^'(5'A2Ar)*l(aA2Ar))) - ^2n)'^\S' 1^2n))) 
X ({{^\s^2N)^\a'/^2N))) - {^\a'/^2N)^\S/^2N))))] } - 

Here, the notation ((...)) stands for the Gaussian average 

^ ^]( )e^^" ^ J dcrtiCT'*(crA2iv)[-R'2Jv(crA2JV,<T'A2Jv)-Q2Jv(o'A2iV,':r'A2iv)]'I'(o''A2iv) 

For the averages {{■■■)) the Wick theorem holds [see Appendix for details], so that 

((¥"((TiA27v)*'^(tT2A27v))) = -Q^^(f72A27V,0-iA2Ar), 
((¥"(criA2Ar)*'5(cr2A2Ar)¥''(a3A2Ar)«''^(CT4A2Ar))) = 

— Q2%{iy i^2N , i^1^2n)Q2n{'^2^2N , Cr3^2N)- 



(3.40) 



(3.41) 



(3.42a) 



(3.42b) 



Inserting Eqs. ( |3.42 ) into Eq. ( 3.40 ) and taking into account Eqs. ( 3.35| ) and ( 3.36| ), one derives after some algebra: 

{p{S)p{S')) = 6{S - S')A2NK2h{SA2N, SA2n) + Ajj^ \k^j,{SA2N, SA2N)k2h{S'A2N, S'A2n) 



k^lj{SA2N, S'A2N)kll{SA2N: S'A2n) + kl%{SA2N. S' A2N)kl],{S A2N : S' A2N) 



(3.43) 



In the large-A^ limit, the kernel K2n{SA2n, S' A2n) 
taken in the b ulk . of . the spectrum is known to obey the 
universal lawsOLjO 



A2Nk2^{SA2N, S' A2n) — A2Nk2N{SA2N, S'A2 

AI^K^],{SA2n,S'A2n) - K\s), 



N 



k2N{SA2N, S'A2n) 



(3.44a) 
(3.44b) 

dtK{t), (3.44c) 



provided the confinement potential is strongo. Here, 
s = and 



sin(7rs) 



ITS 



(3.45) 



Taking into account Eqs. ( |3.44 ) and (3.45), one con- 
cludes that Eq. (3.43) [which is valid for arbitrary s] is 
equivalent to 
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{piS)piS')) = l-K'is) 



2, , dK{s) '■+°° 



ds 



dtK{t) + S{s), 
(3.46) 



that, in turn, is identical to Eq. (1.1a) as one could al- 
ready expect. 



IV. CONCLUSIONS 

The computation—presented above has been inspired 
by the recent papeio by KameiaeiL-and Mezard, followed 
by the further developmentsE^EZlllj, where the problem 
of the fermionic replica rncthod has been revisited. The 
main outcome of Refs.oEjO consists of assertion that 
the replica method in its fermionic u-model elaboration, 
which relies on a selected set of causal saddle pointsEj, 
is capable of reproducing the nonperturbative results 
for the spectral statistics, both in Gaussian invariant 
ensembleaijO in the thermodynamic limit, and in d- 
dimensional disordered conductors with finite Thoulcss 
conductance, g ^ 1, albeit in the large-s domain, s ^ 1. 

In this paper we examined the replica method from a 
different point of view, considering the fermionic and the 
bosonic replicas on an equal footing. The crucial differ- 
ence between the previous studies and the present con- 
sideration resides in using an alternative [to the a-model 
approach] technique to average the replicated partition 
function over the ensemble of random matrices for all 
three Dyson symmetry classes. This technique generates 
such representations for the spectral observables which 



are especially suitable for taking the replica limit n ±0 
in a straightforward way. Performing exact [large-A'^] 
transformations at each stage, we have derived the non- 
perturbative random-matrix-theory fluctuation formu- 
las for the spectral correlation functions in arbitrary spec- 
trum range, including the effect of self-correlation of the 
energy le vels sh owing up in the ^ -functional contribution 
in Eqs. ( |3.12| ), (|3.3l|) and ( ^.46[ ). Let us notice that the 
(5-functional contribution is notpOp academic issue: It is 
this (5-function whose spreadingEZI [in the presence of an 
external perturbation] determines such important char- 
acteristics-pf disordered systemj-a^ distributions of level 
velocitiesEa and level curvaturesEj. 

Reproducing of exact fluctuation formulas for invari- 
ant random matrix ensembles at /3 = 1,2 and 4 through 
the replica method comprises the main outcome of our 
study. It explicitly demonstrates that, at least in the 
context of invariant random matrix models, the replica 
method itself is free from internal subtleties and thus 
completely legitimate, provided one is able to control the 
analytic properties of the replicated partition function 
Eq. ( |1.5| ) as a function of the replica parameter n af- 
ter the averaging over realizations of the random matrix 
Hamiltonian H. The availability of such a control in the 
tj-model-formalizations of the replica method is still an 
open issue. 
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APPENDIX: EQUATIONS (3.27) AND (3.42) 



(i) Consider the average in the l.h.s. of Eq. (3.27a), in which, by definition, 

2N 

W"(a;) = (-l)"-i^Xfc^iri'Hx), 

k=l 

2N 
k=l 

[a,P — 1, 2], see Eqs. ( 3.15| ) and ( 3.1(j ). To compute the required average given by Eq. ( 3.26 ), we write 



2N 



((*"(aiAAr)vI/''(a2A^))) = (-1)"-^ ^ (x^Xfc)^^'^^! A^)^^^ (^2A^), 



(A.l) 
(A.2) 

(A.3) 



where 



{x*iXk) = 



(A.4) 



and M™ is specified by Eq. (3.21). Since the Grassmann integral Eq. (A.4) equals (M ^)kj = Pkj, one derives from 
Eq. (P): 
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27V 



(A.5) 



This completes the proof of Eq. (3.27a). 



(u) Consider the average in the l.h.s. of Eq. ( 3.271: ). Making use of Eqs. ( A.l ) and (|A.2| ), we rewrite the average 



as 



27V 



7',p.q.s — l 

X ^("_-^ViA7v)^ii?\a2Ajv)¥'il^'^(a3Ajv)#f '(a4AA.), (A.6) 



where 



{XrXpX*„Xs) = 



/(n 7n—l 

r. 



m j XrXpXqXs^ 



E2N ^ , , 



/(n m— 1 



f-lprf^sq f-^srf^pq • 



(A.7) 



Inserting this result into Eq. ( |A.6| ) and appealing to Eq. (3.19), we arrive at the r.h.s. of Eq. ( 3.27b| ). This completes 

the proof. 

(iii) Consider the average in the l.h.s. of Eq. (3.42a) in notations of Eq. (3.41). By definition, 

27V 

r(x) = (-i)"-i^x;:[f:'->.-i](^), (A.8) 

2N 

= (A.9) 



fc=l 



a, P — 1, 2], see Eq. ( p.33| ). The required average is 



2A7 



((*"(aiA27v)*''(a2A27v))) = (-l)""' ^ (x;Xfe)[e'">,-i](fTi A27v)[e''-Vfc-i](^2A27v), 



(A.IO) 



where {x*jXk) is given by Eq. (A^) with Mpq defined by Eq. ( 3.37 ). We then conclude that {x*jXk) equals (M ^)kj 



so that 



27V 



((vI/"(aiA27v)^^(a2A2jv))) = (-1)"^' [e'^>,-l](^lA27v)Aifej[e''"Vfe-l](^2A2Ar) = -Q^^((72A27V,aiA27v). 



(A.ll) 



Here, we have used Eg. ( | 3 .36). This completes the proof of Eq. (3.42a). 

The proof of Eq. (|3.42bO " is straightforward and goes along the lines of (ii) and (iii) 
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